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The ground-state properties of spin polarized hydrogen HJ are obtained by means of diffusion 
Monte Carlo calculations. Using the most accurate to date ab initio HJ-HJ interatomic potential 
we have studied its gas phase, from the very dilute regime until densities above its freezing point. 
At very small densities, the equation of state of the gas is very well described in terms of the gas 
parameter pa 3 , with a the s-wave scattering length. The solid phase has also been studied up to 
high pressures. The gas-solid phase transition occurs at a pressure of 173 bar, a much higher value 
than suggested by previous approximate descriptions. 
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I. INTRODUCTION 



The suggestion of Stwalley and Nosanow 1 in 1976 that 
electron-spin-polarized gases of hydrogen (HJ,) could be 
an ideal candidate for achieving a Bose-Einstein conden- 
sate (BEC) state opened an intense experimental search 
which finally led to its first observation in Rb, Na and 
Li in 1995*2 The success attained in cold alkali gases was 
made possible due to techniques like evaporative cooling 
developed previously for confining HJ.. The realization 
of a BEC state of hydrogen was finally achieved in 1998 
by Fried et aL— after overcoming arduous problems like 
recombination on the walls, by working with a wall-free 
confinement, and low evaporation rates by using spin res- 
onance. Also in 1998, Safonov et alX observed a quasi- 
condensate in two-dimensional (2D) HJ, adsorbed on liq- 
uid 4 He. This is one of the best realizations of a 2D 
quantum system since the adsorption energy of HJ. on 
4 He is only ~ 1 K and the adsorbed gas floats approxi- 
mately 8 A apart from the liquid. Further experimental 
work is still necessary to observe in this system signals of 
the Berezinskii-Kosterlitz-Thouless transition, which has 
been recently observed in a trapped gas of Rb confined 
in such a way that atoms can move only within a plane. 5 

Hydrogen is the simplest element and its main prop- 
erties are well known theoretically, starting from the in- 
teratomic interaction which can be computed almost ex- 
actly^ This is significantly different from alkali gases in 
which the interaction is much more involved and, in gen- 
eral, less well known. Its s-wave scattering length a is 
appreciably smaller than the typical values for alkalis, a 
feature that retards evaporative cooling and produces a 
higher transition temperature (50 /iK). Spin-polarized 
hydrogen atoms interact via the triplet potential b 3 E+ 
determined in an essentially exact way by Kolos and Wol- 
niewicz^ and recently extended to larger interparticle 
distances by Jamieson et al.Z The HJ.-HJ, interaction is 



highly repulsive at short distances and presents a shallow 
minimum of ~ 6K at r ~ 4A. The combination of this ex- 
tremely weak attraction and its light mass explains why 
HJ. remains in the gas phase even in the limit of zero tem- 
perature. A measure of the quantum character of a given 
system can be drawn through the quantum parameter j£ 



with e and a the well depth and core of the interaction, 
respectively. According to this definition, 77 = 0.5 for HJ, 
which is the highest value for rj among all the quantum 
fluids (for instance, r\ — 0.2 for 4 He). 

From a theoretical viewpoint, bulk HJ. was first stud- 
ied by Etters et al£ using the variational Monte Carlo 
(VMC) method and by Miller and Nosanow 8 using in- 
tegral equations for computing the multidimensional in- 
tegrals of the variational approach. Both calculations 
showed that the ground state of HJ, is unbound at any 
pressure. Lantto and Nieminen 1 ^ confirmed this result 
using the Euler-Lagrange hypernetted chain equation 
and estimated values for the condensate fraction of the 
gas at different densities for the first time. More recently, 
Entel and Anlaufii carried out a new VMC calculation 
of properties of the gas phase such as the energy, conden- 
sate fraction and excitation spectrum. The heavier iso- 
topes spin-polarized deuterium DJ, and tritium TJ, have 
also been studied theoretically. DJ, atoms obey Fermi 
statistics and their three versions, involving one (Dji), 
two (DJ.2) and three (DJ, 3 ) equally occupied nuclear spin 
states, were analyzed by Panoff and Clarki^ and Flynn 
et al*£ using both VMC and Fermi-hypernetted chain 
theory (FHNC). From the results obtained, they con- 
cluded that DJ.2 and DJ.3 are both liquids at zero pres- 
sure. Their conclusion has been confirmed by Skjetne and 
0stgaarcU 4 -, using a lowest-order constrained variational 
method. On the other hand, microscopic properties of 
bosonic tritium TJ, clusters have been recently studied by 
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Blume et al^ using the diffusion Monte Carlo (DMC) 
method, and their results suggest the use of T| as a new 
BEC gas with the advantage of a nearly exact knowledge 
of the interatomic potential. 

In the present work, we present a DMC study of the 
gas and solid phases of spin-polarized hydrogen. Using 
recent updates of the ab initio Hj-HJ, interatomic poten- 
tial and relying on the accuracy of the DMC method, 
we report accurate microscopic results for energetic and 
structural properties of the bulk system. In the very low 
density regime, the energy is well reproduced by the well- 
know analytical expansion in terms of the gas parameter 
pa 3 , 16 ! 17 ' 18 with a the s-wave scattering length obtained 
by solving the two-body Schrodinger equation with the 
chosen interatomic potential. A relevant result of our 
work is an accurate computation of the gas-solid phase 
transition point which is predicted to occur at pressures 
significantly higher than previous predictions based on 
quantum theory of corresponding states^ and VMC sim- 
ulations.— 

In Sec. II, we report briefly the DMC method and dis- 
cuss the trial wave functions used for importance sam- 
pling in the gas and solid phases. In Sec. Ill, the results 
of the DMC simulations are reported in several subsec- 
tions. In the first one, we review the H|-H| interatomic 
potentials and compare our results at the variational level 
with some previous estimations. The second and third 
parts are devoted to the microscopic results for the gas 
and solid phases, respectively. In the last one, we study 
the gas-solid phase transition point and report results on 
the freezing and melting densities. Finally, Sec. IV com- 
prises a summary of the work and an account of the main 
conclusions. 



II. METHOD 

The DMC method is nowadays a well-known tool de- 
vised to study quantum fluids and solids at zero tem- 
perature. Its starting point is the Schrodinger equation 
written in imaginary time, 

-n^^- = (H-E T MR,t) , (2) 
with an TV-particle Hamiltonian 

^ = -^E^+E^)- (3) 

1=1 l<j 

In Eq. ([2]), E t is a constant acting as a reference en- 
ergy and R = (ri, . . . ,Vn) is a walker in Monte Carlo 
therminology. 

DMC solves stochastically the Schrodinger equation 
© replacing V(R,t) by = V(R,t)ip(R), with 

ijj(R) a trial wave function used for importance sampling. 



In this way equation ([2]) becomes 

- = -DVl${R,t)+DV R {F{R)${R,t)) 

+ (E L (R)-E) <f>(R,t) , (4) 

where D = h 2 /(2m), E L (R) = ipiR)' 1 Hip{R) is the 
local energy, and 

F{R)=2^{R)- 1 V R ^{R) (5) 

is the drift force which guides the diffusion process. In 
Eq. ([¥]), when t — * oo only the lowest energy eigenfunc- 
tion, not orthogonal to ip(R), survives and then the sam- 
pling of the ground state is effectively achieved. Apart 
from statistical uncertainties, the energy of a A-body 
bosonic system is exactly calculated. 

The trial wave function used for the simulation of the 
gas phase is of Jastrow form, ipj(R) — Yii<j f( r ij)> with 
a two-body correlation function f(r) of the form 

f(r) = exp[-6 x exp(-6 2 r)] , (6) 

where b\ and 62 are variational parameters. This form 
has been taken from the VMC work of Etters et al., 9 
who used a Morse potential fitted to reproduce Kolos and 
Wolniewicz ab initio data£ It corresponds to the WKB 
solution of the two-body Schrodinger equation for small 
interparticle distances when the potential is of Morse 
type. 

Simulations of the crystalline bec, fee and hep phases 
have been also carried out; in this case, we use a 
Nosanow- Jastrow model 

N 

^j(R)=MR)H9(ni) , (7) 

i 

where g(r) = exp(— ar 2 /2) is a localizing function which 
links every particle i to a point rj of the lattice. The 
parameter a is optimized variationally. 

The variational parameters b\, 62, and a (|6I7[) have 
been obtained at different densities by optimizing the 
variational energy calculated with the VMC method. For 
example, in the gas phase and at a density p = 0.0079 
A -3 the values are 61 = 82 and 62 = 1-32 A -1 ; b\ in- 
creases with density whereas 62 remains practically con- 
stant. In the solid phase the most relevant parameter is 
a, which increases from a value 0.3 A -2 at melting den- 
sity up to 1.3 A~ 2 at the highest density here studied; 
the Jastrow parameters are kept fixed in all the solid den- 
sity range, b\ — 70 and 62 = 1-32 A -1 . The statistical 
errors of the variational energies are similar to those of 
the DMC results (see Tables [D and HH . 

We use the DMC method accurate to second order in 
the time step At^ and then larger At values than in lin- 
ear DMC algorithms can be used. We have studied both 
the time-step and the mean walker population in order 
to eliminate any bias coming from them. Finally, we 
have analyzed carefully the size dependence of our sim- 
ulations. The calculations on the gas phase have been 
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FIG. 1: HJ.-HJ. interatomic potentials. 
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FIG. 2: Influence of the interatomic potential on the equation 
of state of gas Hj. 



carried out with 128 atoms and some checks with 150 
and 170 atoms have also been made. Using standard tail 
corrections, which assume a uniform system (g(r) = 1) 
beyond r > L/2, with L the length of the simulation box, 
the size dependence of the energy remains smaller than 
the typical size of the statistical errors. The size effects 
are larger in the calculations of the solid phase. In this 
case, we have used 128, 108, and 180 atoms for the bcc, 
fee, and hep lattices. To consider tail corrections in the 
same form as in the gas phase is a rough approximation 
due to the periodic order of the solid. In order to over- 
come this difficulty, we have studied the size dependence 
of the energy at the VMC level where larger number of 
particles can be used. From the VMC results one extracts 
the tail corrections for a given number of atoms and then 
these are added to the DMC energies. It has been veri- 
fied^ that this procedure is able to reproduce accurately 
the experimental equation of state of solid 4 He. 



III. RESULTS 
A. Interatomic potential 

Spin-polarized hydrogen atoms interact via the triplet 
potential b 3 S^, calculated with high precision by Kolos 
and Wolniewicz (KW) in 1965^ Due to the simplicity of 
the H atom it is possible to calculate this potential in an 
essentially exact way. More recently, it has been recal- 
culated up to larger interatomic distances by Jamieson, 
Dalgarno and Wolniewicz (JDW).— The differences be- 
tween the KW and JDW potentials, in the range where 
they can be compared, are rather small, as shown in Fig. 
1. The addition of mass-dependent adiabatic corrections, 
which have been calculated by Kolos and Rychlewski, 22 
to JDW potential, can not be discerned in Fig. 1. 

In the past, only the KW potential has been used in 
the study of the HI gas. Usually, an analytic form was 
assumed and then the free parameters of the model were 
fitted to reproduce the KW data. In this way, Etters et 
al£ used a Morse potential whereas Silvera and Gold- 
mar* 2 ^ proposed a form which is similar to the ones used 



for He-He potentials. The results of these models are 
also plotted in Fig. [1] In the present work, we have used 
the JDW interaction and a cubic spline to interpolate 
between the reported points. The resulting potential is 
plotted as a solid line in Fig. [Tj The JDW data are 
smoothly connected with the long-range behavior of the 
H-H potential as calculated by Yan et alM- The JDW 
potential has a core diameter of 3.67 A and a minimum 
e = —6.49 K (slightly deeper than KW) at a distance 
r m = 4.14A. 

The influence of the potential on the energy of the gas 
at small densities is shown in Fig. O We have carried out 
VMC calculations using the trial wave function © intro- 
duced in the previous Section and the KW and JDW po- 
tentials. The JDW energies are below the KW ones in all 
the density range, reflecting the slightly deeper well of the 
JDW potential. If we use the Morse potential adjusted 
by Etters et al£ to the KW data, the VMC energies are 
significantly worse with respect to the KW ones. This 
manifests the difficulties on fitting a functional form to 
the ab initio KW data; the Morse potential is a bit more 
repulsive than KW and therefore the energies are higher. 
As a matter of comparison, we also show in Fig. [2] re- 
sults from previous calculations. The results from Etters 
et a/3 using the Morse potential are in nice agreement 
with our VMC results using the same potential. The vari- 
ational results of Miller and Nosanow 8 used the KW data 
and arc in close agreement with our present VMC results 
with the same interaction. Finally, results of Lantto and 
Niemineni^ are also reported; they used the KW poten- 
tial and performed an Euler-Lagrange-HNC calculation. 
Their results, restricted to very low densities, are slightly 
better than ours due to their use of an optimized Jastrow 
factor. All these variational results, both ours and from 
previous works, are compared in the same Figure with 
present DMC results with the JDW interaction. As ex- 
pected, the DMC results are below the VMC ones in all 
the density range with a difference that increases with p, 
a predictable feature attending to the fact that the Jas- 
trow factor © corresponds to an analytical form which 
approximates the wave function solution of the two-body 
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FIG. 3: Equation of state of gas H| as a function of the 
gas parameter pa 3 in units of ft 2 / (2ma 2 ) (solid points). The 
triangles correspond to a HS gas 25 and the line to Eq.©. 



problem. 



B. Gas phase 

The gas phase of spin-polarized H at very low densi- 
ties is of special relevance for the field of Bose-condensed 
gases at low temperatures. The case of H is even more 
appealing than alkali gases from a theoretical perspective 
because the interatomic interaction between HJ, atoms is 
very well known, as we have commented in the previous 
subsection. At sufficiently low densities, the equation of 
state becomes universal when it is written in terms of 
the gas parameter x = pa 3 , with a the s-wave scattering 
length. The equation of state of a bosonic gas at low 
densities is given by 



I 1 =J ' T '' 



1 



128 
15V^ 



,1/2 



(8) 



where the first term is the mean-field result^ and the 
second is the Lee-Huang- Yang correction^ the energy 
per particle is written in units of Ti 2 / (2ma 2 ) . 

In Fig. [31 the energy per particle of gas HJ, is com- 
pared to the universal equation of phase ([5]) and to DMC 
results for a hard-sphere (HS) gas from Ref. In order 
to carry out this comparison we have calculated the s- 
wave scattering length of the JDW potential used in the 
present work. The value obtained, a = 0.697 A agrees 
with previous determinations . 15 ' 26 As one can see in the 
figure, the equation of state of gas HJ. coincides with both 
the analytic law © and the HS results up to x ~ 1CP 4 
in agreement with the range of universality determined 
in Ref. Hfl Beyond this value, the energies of bulk HJ. 
clearly separate from Eq. (J5J) , increasing with x faster 
than for HS gas. 

The condensate fraction no, i.e., the fraction of parti- 
cles occupying the zero momentum state, presents also a 
universal behavior in terms of the gas parameter x at very 
low densities. According to the Bogoliubov formula^ 6 - 



FIG. 4: Condensate fraction of gas HJ as a function of the 
gas parameter pa 3 (solid points). The triangles correspond to 
a HS gas 25 and the line to the Bogoliubov approximation ([9} . 
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TABLE I: Results for gas HJ at different densities p: energy 
per particle (E/N), kinetic energy per particle (T/N), pres- 
sure (P), and speed of sound (c). Figures in parenthesis are 
the statistical errors. 



The condensate fraction has been obtained from the 
long-range behavior of the one-body density matrix p(r), 
uq — linv^oo p(r). We have verified by increasing the 
number of particles of the simulation at different densi- 
ties that the size dependence of no is smaller than its 
statistical error. In Fig. |4l we compare our results for 
the condensate fraction of gas HJ. with Eq. ((9]). DMC 
results 25 for no in the HS system are also plotted in the 
same figure. As one can see, the three results are co- 
incident up to x ~ 10" 4 , the same value observed for 
the energy in Fig. [3] Both HS and HJ. show a faster 
decrease with x than the law ©, the departure from it 
being significantly larger for hydrogen, in agreement with 
the same feature observed in Fig. [3] for the energy. 

Spin-polarized H in its gas phase has been studied in 
all the density range up to densities above crystallization. 
DMC results for the total and kinetic energy per particle 
at different densities are reported in Table U In order to 
remove any residual bias from the trial wave function, ki- 
netic energies are calculated as differences between total 
energies and pure estimations of potential energies. The 
energies are positive everywhere, proving its gaseous na- 
ture, and dominated by the kinetic part which is larger 
than the potential energy (in absolute values) at any den- 
sity. The potential energy per particle is negative up to a 
density p = 0.015 A~ 3 , presenting a minimum value of -9 
K at a density p = 0.01 A -3 , and then becomes positive. 

In Fig. we plot the present DMC results for the 



FIG. 5: Energy per particle of gas H| (solid circles) as a 
function of the density p. The solid line corresponds to the 
fit to the DMC energies using Eq. (110 j) . The error bars of the 
DMC energies are smaller than the size of the symbols. 



FIG. 7: Two-body radial distribution functions of the gas 
phase. From bottom to top in the height of the main peak, 
the results correspond to densities 0.002, 0.0067, 0.01, 0.0125, 
and 0.015 A" 3 . 
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FIG. 6: Pressure and speed of sound of gas Hj as a function of 
the density. Left (right) scale corresponds to pressure (speed 
of sound). 



equation of state of the gas. Our results are well param- 
eterized by a polynomial form (e = E/N) 



e(p) = eip + e 2 p 2 + esp 3 + e^p 4 



(10) 



shown as a solid line on top of the DMC results in Fig. 
[5j The best set of parameters is: e\ — 217.0(1.9) KA, 
e 2 = 7.76(9) x 10 4 KA 2 , e 3 = 8.23(12) x 10 6 KA 3 , and 
e3 = 5.1(5) x 10 7 KA 4 , the figures in parenthesis being 
the statistical uncertainties. 

Using the equation of state (fl~0"|) . the pressure is easily 
derived from its thermodynamic definition 



P(p)=p 2 (de/dp) 



(11) 



and from it, the corresponding speed of sound as a func- 
tion of the density 



c 2 { P ) 



1 

m 



dP 
~dp~ 



(12) 



Results for the pressure P and the speed of sound c for 
some values of the density, where specific DMC simula- 
tions have been carried out, are reported in Table [IJ The 
functions P{p) and c(p), derived respectively from Eqs. 
(fTTjl and (fl~2|) . are shown in Fig. [6] 

DMC simulations allow also for exact estimations of 
other relevant magnitudes such as the two-body radial 
distribution function g(r) and its Fourier transform, the 




FIG. 8: Static structure function of the gas phase. From 
bottom to top in the height of the main peak, the results 
correspond to densities 0.002, 0.0067, 0.01, 0.0125, and 0.015 

A" 3 . 



static structure function S(k). With the use of pure esti- 
mators^ 7 - it is possible to eliminate the bias coming from 
the trial wave function and arrive to exact results for both 
functions. The evolution of g(r) with density for the gas 
Hi is shown in Fig. [71 At the smallest density reported, 
g(r) is a monotonic function with the corresponding hole 
consequence of the repulsive core of the interatomic inter- 
action. When p increases g(r) gains structure, with the 
main peak that shifts to shorter distances and increases 
its strength. 

In Fig. results of S(k) at the same densities as 
in Fig. [Jj are reported. The results show the expected 
behavior: when p increases, the strength of the main 
peak increases and moves to higher momenta in a mono- 
tonic way. At low momenta, the slope of S(k) de- 
creases with the density, following the limiting behavior 
limfc— >o S(k) = hk/(2mc) driven by the speed of sound 
c. As expected, the DMC data start at a finite k value 
inversely proportional to the box size L. 

To end this subsection, we show in Fig. [9] the den- 
sity dependence of the condensate fraction from the very 
dilute regime up to densities corresponding to freezing. 
The full set of data is well reproduced using the func- 
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FIG. 9: Condensate fraction of spin-polarized H in the gas 
phase. The line corresponds to a fit to the DMC data us- 
ing Eg . (|13fl . The error bars are smaller than the size of the 
symbols. 



tional form 

noOO - l--^(pa 3 )^-b 1 pa 3 -b 2 (pa 3 ) 2 -b 3 (pa 3 )^ , 

(13) 

which is also plotted in the figure as a solid line on top 
of the DMC results. The values of the parameters in 
Eq. (T3J) are &i = 504(5), b 2 = -1.254(49) x 10 5 , b 3 = 
8.54(55) x 10 5 , and a is the scattering length. 



C. Solid phase 
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TABLE II: Results for solid HJ, at different densities p: en- 
ergy per particle (E/N), kinetic energy per particle (T/N), 
pressure (P), and speed of sound (c). Figures in parenthesis 
are the statistical errors. 
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FIG. 10: Energy per particle of solid HJ, (solid circles) as a 
function of the density p. The solid line corresponds to the 
fit to the DMC energies using Eq. (114[) . The error bars of the 
DMC energies are smaller than the size of the symbols. 



The solid phase of spin-polarized H has been studied 
by using as importance sampling trial wave function a 
Nosanow-Jastrow model 0. The geometry of the lat- 
tice is defined by a proper selection of the lattice sites 
rj around which the atoms are organized according to 
a commensurate solid. There is no experimental mea- 
surement on solid HJ. at the pressures in which we are 
interested in and nothing is known about the form of 
its solid lattice at low temperatures. We have carried 
out calculations of the solid phase at some densities us- 
ing the fee, hep, and bec lattices. Near the melting 
density the bee phase is slightly better and, at higher 
densities, the differences between them are not distin- 
guishable within the statistical noise: at p — 0.0125 
A" 3 , E/N = 33.12(4), 34.30(5), and 33.02(8) K, and at 
p = 0.018 A" 3 , E/N = 76.2(2), 76.2(2), and 76.1(2) K 
for fee, hep, and bec respectively. Therefore, we decided 
to study the solid HJ, properties assuming a bec phase. 
It is worth noticing that the same lattice was used in the 
past by Pierleoni et al£& in the study of solid H at very 
high pressure. 

Some selected results for the total and kinetic energies 
per particle are reported in Table ITT1 The behavior of the 
partial energies, potential and kinetic, is very similar to 
the one obtained for the gas: the kinetic energy domi- 
nates in all the density regime and the potential energy 
is negative at the lower densities and becomes positive 
for densities p > 0.02 A -3 . The full set of results for the 
energy of the solid phase is shown in Fig. |TD] The solid 



line on top of the DMC results correspond to a numerical 
fit obtained using the function 



e (p) = S1P + S3P 3 



(14) 



The optimal values in Eq. fl4|) are si = 1147(6) KA and 
s 3 = 9.57(2) x 10 6 KA 3 . 

From the functional form (fl4|) . and using the corre- 
sponding thermodynamic expressions for the pressure 
(fTTj) and the speed of sound (p~2|) . one can easily derive 
the dependence of these magnitudes on the density. The 
results for both functions are plotted in Fig. [TTJ 

The spatial order of the solid is reflected in the shape 
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FIG. 11: Pressure and speed of sound of solid HJ, as a func- 
tion of the density. Left (right) scale corresponds to pressure 
(speed of sound). 
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FIG. 12: Two-body radial distribution lunctions of the solid 
phase. From bottom to top in the height of the main peak, 
the results correspond to densities 0.0125, 0.015, 0.02, and 
0.025 A -3 . 



of the two-body radial distribution function g{r). Results 
for g(r) in the solid phase for different densities are shown 
in Fig. At densities p = 0.0125 and 0.015 A~ 3 , we can 
compare results for g(r) in the gas and solid phases. As 
one can see, the peaks of the solid are slightly shifted to 
larger distances than in the gas and more importantly, 
and as expected, the secondary peaks of the solid are 
more pronounced. When p increases, the height of the 
peaks increases and moves to shorter distances, like in 
the gas phase. 



50 60 , 70 



l/p(A) 



FIG. 13: Maxwell construction based on plotting the energy 
per particle, E/N as a function of 1/p. The densities at which 
the first-order transition occurs are identified by finding the 
common tangent (solid line) to both the solid (dotted line) 
and gas curve (dot-dashed line). The inset shows the con- 
struction in a wider range of 1/p. 




D. Gas-solid phase transition 

A relevant prediction that can be drawn from the 
present DMC results on the energies of the gas and solid 
phases of HJ, is the location of the gas-solid phase tran- 
sition. In order to determine the transition point and 
the corresponding freezing (pf) and melting (p m ) den- 
sities, we have performed the Maxwell double-tangent 
construction as shown in Fig. 1131 From the common 
tangent to both phases we obtain pf = 0.01328 A~ 3 and 
p m = 0.01379 A~ 3 , which corresponds to a common pres- 
sure at the transition of P = 173(15) bar. The melting 
pressure has proven to be quite independent of the lattice 
used in the simulation since using fee and hep we obtain 
175 and 176 bar, respectively. 

At the transition the kinetic energy per particle of the 
system is discontinuous. The size of this discontinuity 
in other quantum fluids such as He and Ne has received 
the interest in the past from both the theoretical and 
experimental sidesi 29 i 30 ' 31 In the present system we are 
able to accurately measure this discontinuity: in the gas 
side at freezing T/N = 44.0(5) K and in the solid side at 
melting T/N = 51.5(6) K. Therefore, the discontinuity 
amounts to 7.5 K approximately. On the other hand, 
when the system crystallizes the condensate fraction of 
the gas is small but not zero, tiq = 0.04. 

The Lindemann's ratio, defined as 7 = 
yj ((r — T*/) 2 )/ci£, where is the lattice constant, 
can also be obtained from the DMC simulations. At the 



FIG. 14: Two-body radial distribution function at the gas- 
solid phase transition. The solid line corresponds to the solid 
at pm and the dashed line to the gas at pf. 



melting point it is 7 — 0.25, a nearly identical value to 
the one of solid 4 He (7 = 0.26). 

The spatial structure of both phases at the transition 
point is rather different in spite of the small difference be- 
tween pf and p m . In Fig. 1 141 we show results of g(r) for 
both phases at the transition point. As one can see, the 
degree of localization is higher for the solid: the strength 
of the main peak is larger and the height of the subse- 
quent peaks decreases more slowly than in the gas. Nev- 
ertheless, the signature of the solid phase manifests more 
clearly in S(k). In Fig. [15J results of S(k) for both phases 
are compared at densities pt (gas) and p m (solid). High 
intensity peaks located at the reciprocal lattice sites are 
a clear signature of the solid order; they are obviously 
absent in the S(k) of the gas. 



IV. CONCLUSIONS 

The ground-state properties of spin-polarized hydro- 
gen H| have been accurately determined using the DMC 
method. The combination of the accuracy provided by 
the DMC and the precise knowledge of the H|-H| inter- 
atomic potential has allowed for a nearly exact determi- 
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FIG. 15: Static structure factor at the gas-solid phase transi- 
tion. The results correspond to the gas at pi and to the solid 
at p m . 

nation of the main properties of the system, in both the 
gas and solid phases. 

The light mass of the hydrogen atom and the shallow 
well of its interaction force H| to remain a gas in the 
limit of zero temperature. In the very dilute limit, the 
equation of state of the gas is well described by the gen- 
eral expression of a weakly interacting Bose gas. This 
analytical behavior depends only on the gas parameter 
pa 3 , with a the corresponding s-wave scattering length. 
Compared with other Bose-condensed gases, like the al- 
kalis, hydrogen presents the appealing circumstance of 
the accurate knowledge we have of its interatomic inter- 
action. This allows for the use of the real interaction in 
all the density regime, and thus accurate calculations at 



much higher densities are possible. 

When the density is high enough the systems freezes. 
We have studied the energetic and structural properties 
of the solid phase. Near melting the bcc phase is slightly 
preferred over the hep and fee ones. However, the en- 
ergy differences between the lattices are very small and 
become indistinguishable at higher densities. From the 
DMC equations of state of the gas and solid phases, we 
have obtained the gas-solid transition point. At zero tem- 
perature, the transition occurs at P = 173(15) bar. This 
value is significantly higher than previous estimations: 
50 bar, obtained by using the quantum theory of corre- 
sponding states** and 81 bar, from a VMC estimation by 
Danilowicz et al.. 19 It is worth noticing that the transi- 
tion point depends dramatically on the accuracy of the 
theoretical method used for its calculation: if instead of 
using the DMC technique an estimation is performed us- 
ing only the VMC method, one obtains P = 113(17) bar, 
a value significantly smaller than the DMC result. 
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